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Abstract. ATR points were introduced by Darmon as a conjectural con- 
struction of algebraic points on certain elliptic curves for which in general 
the Heegner point method is not available. So far the only numerical evi- 
dence, provided by Darmon-Logan and Gartner, concerned curves arising as 
quotients of Shimura curves. In those special cases the ATR points can be 
obtained from the already existing Heegner points, thanks to results of Zhang 
and Darmon-Rotger-Zhao. 

In this paper we compute for the first time an algebraic ATR point on a 
curve which is not uniformizablc by any Shimura curve, thus providing the 
first piece of numerical evidence that Darmon's construction works beyond 
geometric modularity. To this purpose we improve the method proposed by 
Darmon and Logan by removing the requirement that the real quadratic base 
field be norm-cuclidcan, and accelerating the numerical integration of Hilbert 
modular forms. 



1. Introduction 

Let F be a totally real number field and let E/F be an elliptic curve of conductor 
9t. Denote by L(E/F, s) the Hasse-Weil L-series attached to E, which is known to 
converge in the half plane Si(s) > 3/2. Let us assume thorough this note that E 
is modular; that is to say, that L(E/F, s) equals the L-series of a Hilbert modular 
form over F of weight 2 and level 9L Thanks to the modularity theorems of |Wil95j , 
[BCDTOlj and [SWOlj . E is known to be modular if either F = Q, or if [F: Q] > 1 
and it satisfies certain mild conditions on the reduction type at primes above 3. The 
L-series L(E/F, s) admits analytic continuation, and the Birch and Swinncrton- 
Dyer (BSD) Conjecture predicts that its order of vanishing at s = 1, called the 
analytic rank of E/F, equals the rank of the group of F-rational points E(F). 

In this context, BSD Conjecture is known to hold in analytic rank or 1 provided 
that E satisfies the following Jacquet-Langlands hypothesis: 

(JL) Either [F: Q] is odd or there is a prime p in F such that ord p (9t) is odd. 

Theorem 1.1 (Gross-Zagier, Kolyvagin, Zhang). Let E be a modular elliptic curve 
over a totally real number field F satisfying (JL). 7/ord s= i L{E/F 1 s) < 1 then 

ord s= i L(E/F, s) = rank(E(F)). 

In analytic rank 0, BSD is also known for modular elliptic curves not satisfying 
(JL), thanks to the work of Longo |Lon06j . However, in analytic rank 1 (JL) cannot 
be dispensed with at the moment, because the construction of non-torsion points 
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relies on the existence of the so-called Heegner points. Indeed, if E satisfies (JL) 
then it is geometrically modular: there exists a non-constant i^-homomorphism 



from the Jacobian of a suitable Shimura curve X defined over F onto E. Shimura 
curves are endowed with CM points, which are defined over ring class fields of 
quadratic CM extensions K/F. The projection of CM points via tie gives rise 
to Heegner points on E, whose arithmetic behavior is linked to the corresponding 
L-series of E thanks to formulas of Gross-Zagier and Zhang. 

On the other hand, if E does not satisfy (JL) then it is not known to be geo- 
metrically modular unless it is a Q-curve; i.e., a curve isogenous to all of its Galois 
conjugates. As a consequence of Serre's modularity conjecture Q-curves admit Q- 
parametrizations from classical modular Jacobians, and this has been exploited in 
[DRZ12] in order to construct Heegner points and prove BSD in analytic rank 1 
for some Q-curves not satisfying (JL). But BSD in analytic rank 1 seems to remain 
intractable for elliptic curves which are not Q-curves and do not satisfy (JL). In- 
deed, since they are not geometrically modular the Heegner point method sketched 
in the previous paragraph cannot be applied in this setting. 

The Heegner point construction constitutes the only known procedure for sys- 
tematically manufacturing algebraic non-torsion points on elliptic curves. However, 
several conjectural constructions have emerged in the last years under the generic 
name of Stark-Heegner points, or also Darmon points as the first such construction 
was introduced in [DarOlj . Variants of this initial construction applying to several 
different settings have been proposed since then, for instance in |Das05j . [Gre09j . 
|LRV09j . [Gar 1 la] , and |GRZ12j . The leitmotif of these methods is the analytic 
construction of algebraic points on ring class fields of quadratic extensions K/F 
which, unlike the classical not CM. 

This note deals with the effective computation of a type of Darmon points known 
as ATR points, which where introduced in |Dar04| Chapter VIII]. To explain the 
terminology, recall that a number field is said to be almost totally real, or ATR for 
short, if it has exactly one complex non-real archimedean place. 

Let K be a quadratic ATR extension of the totally real field F. For an ideal c of F 
denote by R c C K the order of conductor c, and H c the ring class field corresponding 
to R c . Darmon associates to any optimal embedding <p: R c M2(Of) a point 
P v € E(C), called an ATR point, which is conjectured to be defined over H c . 
Moreover, by analogy with the formulas of Gross-Zagier and Zhang, its trace down 
to K is believed to be non-torsion if and only if ord s= i L(E/K, s) = 1. 

An algorithm for computing ATR points in the particular case where F is real 
quadratic and E has conductor 1 is given in |DL03j . These elliptic curves do not 
satisfy (JL), so that the Heegner point construction is not available in general. Both 
the definition of the points P v and the conjectures of Darmon concerning them will 
be recalled in Section [5J For the moment, it is enough for us to mention that they 
are the image under the Weierstrass map C/A^; — > E(C) of a complex number of 
the form 



where w is a certain differential 2-form on the Hilbcrt modular surface SL2(C_f)\'H 2 . 
The limits of integration depend on the embedding <p, but they are not uniquely 
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determined: for a given tp there are many possible choices for ai,bi,Ci,di. The 
Fourier series of uj is explicitly computable from E, and term by term integration 
of a truncation leads to a numerical approximation to J v . The rate of convergence 
depends essentially on the imaginary parts of the limits, and this turns out to be 
the main computational restriction of this method. 

The algorithm outlined above was used in |DL03j to obtain numerical evidence 
towards Darmon's conjectures. More concretely, ATR points on three concrete 
elliptic curves were computed, and they were checked to be (up to a certain nu- 
merical precision) multiples of the corresponding generators of the Mordeil-Weil 
groups. Calculations of the same kind were performed in [Garllbj for one more 
curve. However, computational limitations restricted them to elliptic curves which 
all happen to be geometrically modular. In this case, the BSD Conjecture im- 
plies that ATR points in these curves should be related to the already existing 
Heegner points. Actually, in the case of Q-curves (as they are all examples consid- 
ered in |DL03j ) a precise relation between Heegner and ATR points is conjectured 



In this note we speed up the algorithm devised by Darmon and Logan by im- 
proving its two main bottlenecks. Namely, the computation of integrals of Hilbert 
modular forms, and the determination of limits dj, Cj, di in Q1.2[) having the high- 
est imaginary part possible. This allows us to gather more numerical evidence in 
support of Darmon's conjectures, by calculating ATR points on elliptic curves which 
were not computationally accessible using the algorithm in |DL03| . In particular, 
we have been able to compute for the first time an ATR point of infinite order on 
a non- geometrically modular elliptic curve. 

More concretely, the contents of the article are as follows. In Section [5] we 
review the definition of the points P v and Darmon's conjectures on their arithmetic 
properties. We also sketch the algorithm of Darmon and Logan for their explicit 
computation. 

In Section [3] we present the algorithm for speeding up the computation of in- 
tegrals of Hilbert modular forms. The idea is to use the fact that the limits are 
invariant under the group SL2(Of) in order to transform the given integral into a 
sum of integrals whose limits are uniformly bounded away from the real axis, the 
bound depending only on F. It is worth remarking that this algorithm does not 
exploit any particular property of the integrals involved in ATR points, and there- 
fore it may be of independent interest for computing integrals of Hilbert modular 
forms in other contexts. 

In Section [4] we comment on a trick that can sometimes accelerate the compu- 
tation of ATR points. The procedure for computing the limits ai,bi,Ci,di in (|1.2j) 
for an embedding tp involves at some point the calculation of a continued fraction 
expansion of an element in F. We exploit the non-uniqueness of continued fractions 
to attach to a given tp limits a%, hi, c,, di with as high imaginary part as possible. 

Finally, in Section [5] we use Darmon-Logan's algorithm together with the im- 
provements of sections [3] and [4] to enlarge the pool of elliptic curves on which Dar- 
mon's conjectures have been numerically tested. Arguably, the most interesting 
among them is the curve E^qq given by the equation 



in [DRZ12I §4.2]. 



y 2 - xy - uy = x 3 + (2 + 2w)x 2 + (162 + 3w)x + (71 + 34w), 
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because it is not a Q-curve. We have computed an ATR point corresponding 
to the field K = F(\/9144a; + 98577), and we have numerically checked that it 
coincides with a multiple of the Mordcll-Weil generator of £509 (if). Since -E509 is 
not geometrically modular, such point docs not seem to be explained by the presence 
of Hccgner points. This gives experimental evidence that Darmon's construction 
leads to algebraic points that arc genuinely new, not attainable by classical methods. 

Finally, it is worth mentioning that ATR points are also the base of an algorithm 
by L. Dcmbele |Dem08j for computing equations of elliptic curves with everywhere 
good reduction attached to Hilbert modular forms of level 1. The authors hope that 
the algorithm presented in this note can be useful for this purpose, and that it may 
lead in the future to a systematic computation of such equations using Dembcle's 
method. 

Acknowledgments. The authors are thankful to Henri Darmon and Victor Rotger 
for initially suggesting the problem and for many helpful conversations, and to the 
anonymous referee for many valuable comments and suggestions. They are grate- 
ful to the Max Planck Institute for Mathematics for the hospitality and financial 
support, and for making available its computational resources, crucially needed for 
part of this note. Part of this work was also carried out in the facilities of Centro 
de Ciencias de Bcnasque Pedro Pascual during the Summer of 2011. The authors 
received financial support from DGICYT Grant MTM2009-13060-C02-01 and from 
2009 SGR 1220. 

2. Computation of ATR points 

Let F be a real quadratic number field of discriminant D and narrow class 
number 1. Write Of for its ring of integers and set T = SL2(C_f). Wc denote by 
Vq,Vi the embeddings of F into R. For an element iff we may write Xi instead 
of Vi(x), and |x| for the norm of x, given by xqX\. Recall that T acts discretely 
on T~L 2 via vq x v\. The analytic variety V \ TL 2 can be compactified by adding one 
cusp, which gives rise to the Hilbert modular surface X attached to T. 

Let K/F be a quadratic ATR extension. For an ideal c of F we denote by H c 
the ring class field corresponding to the order i? c of conductor c of K. Wc assume 
without loss of generality that vq extends to a complex place of K and v\ extends 
to a pair of real places of K. We fix an extension of «o to Q C C, which we use to 
to identify K and its extensions with subficlds of C. Recall that an embedding of 
©F-algebras p: R c ^> M 2 (F) is said to be optimal if (p(K)n M 2 (O f ) = <p(Rc)- We 
denote by £ c the set of such optimal embeddings. 

Let E/F be an elliptic curve of conductor 1. In this section we review Darmon's 
construction, which attaches to each optimal embedding ip a point P v G E(C) that 
conjecturally belongs to E(H C ). There are several equivalent ways of defining P v . 
For instance, a nice geometric definition in terms of a non-algebraic analogue to the 
Abel-Jacobi map is given in |DRZ12| §2.1]. However, for computational purposes 
the original definition of |DarQ4[ Chapter VIII] , or rather the subsequent refinement 
of [DL03] are better suited. Key to the approach of |DL03j is the definition of 
certain semi-definite integrals of Hilbert modular forms, whose existence and main 
properties we will take as a black box. 

2.1. Semi-definite integrals of Hilbert modular forms. Let / e S2 (T) be a 

Hilbert modular form. Recall that / has a Fourier expansion indexed by totally 



COMPUTATION OF ATR DARMON POINTS 



5 



positive elements of Of- Actually, the Fourier coefficient corresponding to n € Op 
only depends on the ideal (n) generated by n, and the expansion is of the form 

(2.1) /(r ,r 1 )= « (ll) e 2 " ( ^ ro+ ^ ri) , V^W, 

where 5; = and Sf is a totally positive generator of the different ideal of F. 

Let us assume from now on that all Fourier coefficients a(„) are rational numbers. 
The reader can refer to |GRZ12| §2.4] for the definition of ATR points when the 
Fourier coefficients generate a number field of degree > 1, in which case they belong 
to some higher dimensional modular abelian variety. 
The differential form 

Uf = — 7=-/(ro,ri)dToari 



D 

is invariant under the action of T and extends to a holomorphic form at the cusp, 
thus defining a holomorphic 2-form on X. The expansion in Equation (|2.1[) is useful 
for computing integrals ofcjf. Indeed, for xq , Xi , yo , y\ £ H we have that 



(2.2) 



In the definition of ATR points the key role is not played by ujf but instead by the 
non-holomorphic differential 

(2ni) 2 ( \ 
= r— (/(To,ri)dr dri + f(u T , uiTi)d(u T )d(urfi) j , 



D 

where it is a fundamental unit in F such that uq > 1 and u\ < — 1. The differential 
form ui^ is also T-invariant, and it follows from the definition that it is invariant 

under the action of the matrix 7„ = ^ ^ . Therefore, if we let 
f = { 7 e GL 2 (O f ) : «o(det(7)) > 0} 

we find that 

rvo fy\ rfyo riyi 

ujf = / ojf, for all 7 G f. 



/ - / / "7 

x Jx 1 Jix Jyx 1 



Here 7 acts on the outer limits (reps, inner limits) through vo (rcsp. vi). 
Let 



A+ = : 5eH 2 (X,Z)^ C 



be the period lattice of uj^ . 

Theorem 2.1 (Darmon-Logan). There exists a unique map 
H x ¥ X (F) x P 1 (F) — ► C/A+ 

( t o,ci,c 2 ) ' — > r° f^tf 

satisfying the following properties: 

(i) r° /;? = r i: v «« 7 e f , 
(a) r r: 4 + r c <4 = r c <4> 



(i 
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For the existence of such map we refer the reader to |DL03| §1]. Uniqueness is 
proved in |DL03| §4], and it follows from repeated application of properties (0), (HU 
and (jiul) . Since this also leads to an algorithm for computing the map, we review 
the proof in the next section. 

2.2. Computation of semi-definite integrals via continued fractions. Given 
bo,bi,..., 6„ G Of the (finite) continued fraction [bo, b\, . . . , b n ] G F is defined, 
inductively, as 

[b ] = b , [b ,h] = 6 + t-, • • • , [bo,h,. ..,b n ] = [b , [h, . . .,b n ]]. 
Let [bo, bi, . . . , bk] = with pk, qu G Of coprime. It is well known that 

(2.3) p k q k -i-Pk-iqk = {-^) k ~ r - 

Since Op is infinite and F has trivial class number, by a result of Cooke- Vaserstein 
[Coo76j every element c E F can be written as a finite continued fraction. Sec 
[GM12] and 2] for an effective version of this result. 

Two cusps c\, C2 G P^-F) are said to be adjacent if ci =7-0 and C2 =7-00 
for some 7 € T. One can join every cusp c € F with co by a sequence of adjacent 
cusps. Indeed, if c = [bo, b\, . . . , b n ] then the sequence 

Po Pi Pn 
00, —,—,...,— = c 

qo qi q n 

has this property thanks to (|2.3p . Using this fact and property Q every integral 
f T " Ici w / can ^ e wr hten as a sum of integrals of the form f T ° f^. °° . Thanks 
to (0) they are of the form J T J °° u>f, and can be computed as follows: 

p — — r — l pT—1 roc 

=1 ' L -W L 

Integrals of the form J °° u>^ , with t\,t% € T-L, can be computed by taking T3 £ H 
and writing 

PT2 POO [>T2 /*T3 t>T2 |*00 /> — /• — /*T2 fOC 

(2.4) / / w ; = / / / .+ = /;/%++/ / 4, 

for which formula (|2.2p can be used. 

2.3. Definition of ATR points. Under the assumption that E is modular there 
exists a Hilbcrt modular newform Je G ^(r) such that L(E/F, s) = L(fE, s). The 
Fourier expansion of Je can be explicitly computed as follows. For a prime ideal 
p C F, let a p := \p\ + 1 - #E(0 K /p), where |p| denotes #0 F /p. For arbitrary 
ideals ncf the integer a n is defined by means of the identity 

nCF p prime 
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and the Fourier expansion of Je is given by 

(2.5) f E (r ,r 1 )= a (n) e 2 " ( ^ ro+ ^ Tl) , To^eH. 

neo% 

Let if: R c M2(F) be an optimal embedding. Since K is ATR the group 

T v = {7 G <p(0 F ) : det( 7 ) = 1} 

has rank 1; let 7^ be one of its generators. Since Vq extends to a complex place of 
K, the action of K x on H by means of vq ° ip has a single fixed point tq. Let J v be 
the quantity in C/Ajr defined as 

/to r'/ipoo 
/ W /' 
J oo 

Let S H°(E, fi 1 ) be a differential which extends to a smooth differential on 
the Neron model of E over Of- Let 



A,= 



lJ z v iM ■ Z e H 1 (Ej , Z) I 



and let (rcsp. Aw ) be a generator of Aj n K (resp. Aj n zR). 
Conjecture 2.2 (Oda). There exists a integer c such that -ttAj C Ao- 

Granting Conjecture 12.21 and denoting by 77: C/Ao — > -Eo(C) the Weierstrass 
parametrization map, the ATR point is defined as 

The group T acts by conjugation on £ c , and the set of equivalence classes £ C /T has 
cardinal [H c : K] (cf. |Uar04l §8.5]). 

Conjecture 2.3 (Darmon). The point P v belongs to E Q (H C ). Moreover, the point 

Pk= £ p v 

<pe£ c /f 

belongs to Eq(K) and it is non-torsion if and only if ord s= i L(E / ' K, s) = 1. 

Note that under the assumptions of this section, namely the conductor of E/F 
is 1 and K/F is ATR, the L-function L(E/K,s) has sign — 1. Thus the condition 
ord s= i L(E/K, s) = 1 is equivalent to L'(E/K,1) 7^ 0, which can be numerically 
checked. 

3. Integration of Hilbert Modular forms 

When evaluating (|2.2j) one can collect the elements n £ Op modulo powers of 
u 2 , because the Fourier coefficient corresponding to n only depends on the ideal 
(n). The sum corresponding to n ■ (u 2 ) is then 



00 

k=—oo ^ 
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For a fixed k, the modulus of each exponential term (after multiplying out the 
brackets) is of the form 

(3.1) e - 2 -Q™o fc +^«o 2fc ) j with r € {3( 2 / ),3(x )},se {3( yi ),3(a;i)}. 

It is easy to see that (|3.1j) has a single maximum, when viewed as a function of k, 
and that the maximum value afforded by the four exponential terms is bounded by 

M n = e vn , 
where e = e(xo, yo, xi, Hi) is defined by: 

(3.2) e(x ,yo, xi,yi) 2 = min{3(xo)3(xi), 9f(xo)3f(j/i), 3(yo)3(yi), S(yo)S(xi)}. 
Moreover, one easily checks that S n is dominated by a geometric series and that 

S n < 7-7^ , for n » 0. 

|n| 

This estimate allows us to know a priori how many terms need to be considered 
in order to obtain a prescribed accuracy. We observe that the speed of conver- 
gence of expression (|2.2j) depends on the limits Xo, yo-, x\, yi through the quantity 
e(xo, yo, x\, y\). The main result of this section is the following. 

Theorem 3.1. There exists a constant e F , which depends only on F, such that 
for every {xo,yo,X\,yi) G "H 4 and for every e < e F the integral f x ° J" 1 ujf can be 
expressed as 

(3.3) / / Uf = / + / / cJ f , 

with e(oi, hi, Ci, di) > eo for all i = 1, . . . , n. 

Observe that in the integrals considered in Theorem 13.11 the four limits of in- 
tegration lie in TL. By (|2.4p integrals of the form fj° uj^ are involved in the 
computation of ATR points. One can choose any y\ G TL whose imaginary part is 
large enough to satisfy that 

rVo roo rya ryi rya roo 

/ / w / = / / w / + / / w /' with s ( a; o)3(yi),3(yo)S(yi) > e F , 

J xq J xi J xq J x± J xq J y± 

and apply Theorem 13.11 to the integral JJ 1 oj^ . 

We devote the rest of the section to prove Theorem 13. 1[ which as we will see can 
be made effective and algorithmic. 

It will be useful for us to regard F as a subset of M 2 by means of vq x v\ . Let 
|| ■ || denote the norm on R 2 given by 

||(x ,xi)|| = max{|x |, 

The basic ingredient in the proof of 13. II is the following classical result. 

Lemma 3.2. There exists a constant C F , only depending on F, such that for each 
x € R 2 and for each < 8 < 1 there are elements c, d G Ok, with c ^ 0, such that 

Hec + dll^*, ||c||<%. 

6 
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Proof. This is |Fre90| Lemma 3.6]. We rewrite the proof in an algorithmic fashion 
in order to give an approximation to Cf (see Remark 13.31 below) . as well as to give 
an algorithm to find the elements c and d. 

Consider a fundamental parallelogram P for Of as a subgroup of R 2 . Let 
Ui, . . . , Un be boxes of side 5 that cover P. It is easy to see that there is a constant 
N', depending only on F, such that N can be taken to be < N' /S 2 . 

For each positive real r, consider the set 

S F (r) = {ceO F | ||c||<r}. 

A well-known result in Ehrhart theory (see e.g. |BR07( Theorem 2.9]) implies 
that there exists a constant Cf > 0, which depends on F but not on S, such that 

#Sf ( Cf) > £ > w . 



v 26 J S 2 

Consider now an ordering {c n } n >i of the elements of Sf(Cf / (26)). For each of the 
c„, find d n € Of such that 

c n x + d n e P, 

and set i(n) to be the integer such that c n x + d n 6 U^ n y By the pigeonhole 
principle the sequence {i(n)} n will have a repetition, say i{n{) = i(ri2). Therefore 

\\(Cni ~ c n2 )x + (d ni - d n2 )\\ < S, 

and 

ii n n n M M Cf Cf 

\\c ni -c„ 2 || < ||c ni || + ||c n2 || < = — ■ 

The seeked elements are thus c = c„_! — c„_ 2 and d = <i n _i — d n —2- D 

Remark 3.3. In our applications, the parameter S will be small enough so that the 
quantity N in the above proof can be approximated by area(P)/<5 2 . Also, |BR071 
Theorem 2.9(b)] gives in this case that 

/ r \ r 2 4 
\2oJ Ao z area(P) 

Therefore, it is enough for #S , _f(^j) to be larger than N that 



,.2 



4 area(P) 



4<5 2 area(F) 5 2 

This yields an approximate value of Cf ~ area(P), which is good enough for 
our purposes. Note that this area is easily calculated: if Of = Z © Zw, then 
area(P) = \wq — w\\. 

We define 

_ u 
eF ~ C F {l + u 2 Y 

In the next lemmas we prove that this can, indeed, be taken as the constant ep of 
Theorem 13.11 

Lemma 3.4. Let z = (zq, z\) be an element inl-L xH. There exists a matrix 7 £ T 
such that 

(3(7oz )3(7iZi)) 1/2 > e F . 
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Proof. Let Zj = r.j + isj and let r = (tq, r\). If sqs\ > e 2 F one may take 7 = 1 and 
there result is obvious, so assume from now on that sqsi < t 2 F . By Lemma 13.21 for 
each < S < 1 there exist d , d' £ Of with d 7^ and such that 

||c'r + d'||<5, ||c'|| <^. 

We have that 



((c r + d ) 2 + d 2 s 2 ) ((dm + d' x ) 2 + c? S ?) < (V + 



£ 2 J \ S 2 

Let g = gcd(c', d!) and let n = Nm^ /q(<?) ■ If we let c = d / g and d = d'/g then 

n2 2 \ / ,-r2 2 



((c r + d ) 2 + c 2 s 2 ) ((cm + rfx) 2 + c 2 s 2 ) < -1 (V + 

<(V + ^#) (^ 2 + ^i 

Since gcd(c, d) = 1 there exists 7 € T having (c, d) as bottom row and 

3(7o zo)Q(7i z i) 



c z + d 1 2 |ciZi + di\ 2 

S Sl 



> 



((c r + d ) 2 + c 2 s 2 ) ((cm + di) 2 + c 2 s 2 ) 

SqSi 



(s 2 



We choose 5 that maximizes this expression. The optimal value for S turns out to 
be 

S= (C 2 F s oSl ) 1/4 . 

Note that S < 1 since soSi < ep. With this value of S the above inequality gives 

ccv \»r ^^ 1 /2 ^ 1 G(s ,si) 

(9(70^0)9(71^1)) > ^-77 v 

2Cf -4(so, si) 

where G(sq, si) and A(so 7 si) are the geometric and arithmetic means, respectively. 
Of course, for this quantity to be not too small we should ensure that the ratio 
G/A is not too small. That is, sq and si should be close. This can be done by the 

action of the matrix 7« ~ - lj j which guarantees that: 

(3.4) u 2 <^<u 2 . 

so 

Therefore one obtains (the worst case is when the ratio is at any extreme): 

(3.5) > _^L_, 
V ' A(s ,si) ~ 1 + u 2 ' 

and the result follows. □ 

Remark 3.5. A similar argument to that of Lemma |3 .41 shows that for any (zq 1 z\) € 
H x T-L there exists 7 € T such that 



(3(7o,o)3(7i^)) 1/2 > r V 
C F (1 + u ) 
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Indeed, in this case one can improve p.4|) to 

-l ^ Vi 
u Q < — <u 

yo 

by using the action of 7 U = ^ jj e f. We will apply this remark to the 

computation of ATR points, because the integrals JJ° are T-invariant, and 

then we can take e F in Theorem 13.11 to be 



(3-6) e F = ,\ — r. 

C F {1 + uo) 

Given x,y g H denote by p(x, y) the geodesic in T-L joining x and y. We also 
let d{x, y) be the hyperbolic distance between x and y. This distance is invariant 
under the action of SL2 (R) , and it is given by the formula 

coshd{x > y) = 1 + ^(x)lly 

Observe that 

(3.7, «,,,)> COS-- (l + (3 gf,'gn )2 

For Xq, J/oj ^lj 2/1 S ^ we denote by ^} } the image under the quotient map 
H 2 U {00} -> T \ (H 2 U {00}) of the region p(x ,y ) x CH X 71. Since 

SL 2 (R) acts by isometries on "H we have that { 1° xi } = { 7^0 ™\ } for all 7 e T 
and the area of { ^° ^ } is given by d(x , y )d(x 1 ,yi). 

For a given eo < e F , let e\ := y/eae F . This auxiliary quantity will be used in the 
proof of Theorem 13. II We also define quantities do, di and d m i n as 



cosh" 1 ( 1 + <y<y€F J, tl \ ^ ) for i g {0, 1}, and d min = min{di, d 2 }. 



Lemma 3.6. Let Xq, yo, Vi ^ % be such that e(xo, yo-, x±, y±) < eo < e F . Then 
fx° Ixi w / can t> e written as a sum of either two or three integrals along sub-regions 
of { x° x\ } with zero-measure intersection, in which the first one is of the form 
Ix° fx' 1 w / w ^ x 0i Vol x 'ii Vi satisfying the following properties: 

(1) e(x' ,y' ,x[,y{) > e 

(2) for each i g {0,1} either = (xi,yi) or d(x-,y-) > d mm ; 
Proof. Actually we will prove that x ,y , x[, y[ can be chosen to satisfy 

(3.8) 3(4)304) > 4, H4My'i) > e?, > 4, WoWO > 4- 

Observe that Q3.8P implies condition (1) of the lemma. 

By Lemma 13.41 and transforming the limits xq, yo, y\ of the integral under an 
appropriate clement of T, we can assume that 

(3.9) $S(xoMxi) > 4. 

If 3(x )3(yi) > el, 3(xi)3(y ) > 4 and 3(y )S(yi) > e\ then e(x 0) 2/o,2i,2/l) > 
eo. So we can suppose that one of these conditions is not satisfied. We divide the 
proof in three cases, that we label (a), (b) and (c), according to the condition which 
is not satisfied. 
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(a) If 9(cco)9(j/i) < e 2 . Then we break the integral as 

ryo ryi ryo rti ryo rVi 

(3.10) / u> f = / u f + / u f , 

where t\ £ p(xi,yi) satisfies that 

(3.11) 9(*i)9f(a;o) = e?. 

e 2 e 2 e 2 

Observe that this is possible because $3(xi) > cj^y > an d ^(z/i) < a(g ) > 

so that the geodesic between xi and j/i in Hi intersects the line $t(z) = a ^ ) , 
and we take ii to be the intersection point. 

If 9(j/o)^(*i) > e o then multiplying this condition by (|3.9p and using p. lip 
we obtain that 3(xi)$3(j/o) > e o an d dHJ) holds. Moreover, dividing (|3.9p by 
(pTTTj) wc find that 

3(xi)/3(ti) > (ejr /ei) 2 
so that, in view of (|3.7p . the first integral satisfies condition fl?). 
If 3(y )9(*i) < e 2 then we break ([3"Hi]) as 

/■l/o /-j/i /-to /-ti /*yo /-I/o /-I/i 

(3.12) / u>f= "f+ / <"/ + / / w /> 

Jxq Jx± J Xq J X\ J to J Xi J Xq J t\ 

where we take £o £ p(%o, Vo) such that 

(3.13) S*(<o)$*(*i) = eg- 

e 2 e 2 e 2 

This is possible because 3(xo) = 5 ^ > 5 ^ and 3(y ) < ■ Now 
multiplying ()3.9p and (|3.13p and using (|3.1ip we obtain that S(xi)3(£o) > Co- 
in addition, dividing (|3.9p by Q3.13P we find that 

9(x )/3(i ) > (ejr /e ) 2 , 

which implies condition (2) for the first integral. 

(b) If 3(xi)3(?/o) < £?• This is analogous to the first case: we break the integral 
as 

ryo fVi rto fVi fVo rVi 

(3.14) / / Uf = / / u f + / / W/) 

where to € p(xo,yo) satisfies that 53(io)5(xi) = e 2 . If 53(fo)53(j/i) < eg then we 
break further the integral as 

ryo ryi rto rti rt ryi ryo ryi 

(3.15) / / u f = / / + / "f+ <*f, 

J Xo J Xl J Xo <J X\ J Xo <J t\ J to J XI 

where we take ii € p(xi, yi) such that S3(fo)3(£i) = eg- 

(c) If !J(2/o)3(?/i) < £?• In this case we can assume that 

(3.16) 3(x )3(yi) > e 2 and 9(xi)3(i/o) > e?, 

because otherwise we are in case (a) or (b). Then we can also suppose that 
(3-17) ^(yoMvi) < e 2 , 

because otherwise e(xo, yo, xi, yi) > eo- Then we can break the integral as 

ryo ryi rto ryi ryo ryi 



'to J Xi 
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where to € p(xo,yo) satisfies that 
(3.18) 3(*o)3(2/i) = eo- 

Observe that 

9(*o)S*(*i) = > > e 2 , 

S(yi) 

so that the first integral satisfies the conditions of the Lemma, because from 
(1316)) and (13T81) if follows that ^(x )/^(t ) > (ei/eo) 2 = (eF/ei) 2 . 

□ 

Proof of Theorem 13.11 Define the quantity e(I) as 

/•2/0 ryi 

e(I) = e(x ,yo,xi,yi), if /= / / w/. 

We can compute an expression / = I± + ■ ■ ■ + I n with e(ij) > eo as in (|3.3I) by 
repeated application of Lemma 13.61 To be more precise, we can apply the following 
algorithmic procedure. 

We initialize a void list V and a list W — [I]. We denote by Wo the first element 
in W. If e(Wo) > eo we append Wo to V and we remove Wo from W. If e(Wo) < eo, 
we apply Lemma 13.61 to write Wq = Iq + J\ + • • • + Jk with e(7o) > eo and k < 2. 
We append Io to V, we remove Wo from and we append Ji, . . . , Jk to W. Then 
we repeat the process to the first element Wo in the list W, until W is void. 

Observe that if this process finishes in a finite number of steps, then at the end 
the list V contains integrals Iq, 1%, ■ ■ ■ ,I n with I = Iq + ■■■-{- I n and e(Ii) > eo for 
alii — 1, . . . , n, as desired. What remains to be shown, then, is that this procedure 
cannot be repeated infinitely many times. As we shall now see, this follows from 
properties (2) and (3) of Lemma 13.61 

First of all, observe that I is the integral of w/ along the region { ^° %\ }. Applying 
Lemma [3.6l to a certain integral Wo — w / amounts to give a decomposition of 

the region { *° *\ } Q { x° x\ } into a union of regions with zero measure intersection. 
This decomposition can be of the following 4 types, which correspond to (13.121) . 
([535]) . (pnU|) and respectively: 

(i) { Z r\} = { ll t\ } U { l\ } U { *° g } with d(r i5 ti) > d rain for * - 0, 1, 

(H) { r°o *[} = { r°o r! } U { *° ^ } U { ? ° ^ } with d( r<1 ti) > d mi „ for S = 0, 1, 

(iii) { °» S r [} = { s r ° *\ } U { s r ° H } with d(r u h) > d min , or 

(iv) { r° S r{} = { % S r \}U{ Z r\ I with d(r ,t ) > d min . 

Each time that an application of Lemma 13.61 gives a decomposition of type (i) 
or (ii), the first term is a subregion of { 1° } of area at least d 2 ^. Since the area 
of { x° x\ } is finite, this cannot happen infinitely many times. Therefore, in the 
algorithmic process described above we can assume that, after a finite number of 
steps, all applications of Lemma 13.61 give rise to decompositions of types (iii) or 
(iv). 

At this stage, the algorithmic procedure applied to an integral Wo is as follows. 
Let Ro = {ro rl} be the region associated to Wo, and let A(i? ) = d(r ,s ) + 
d(ri, Si). As application of Lemma T3.6I one obtains a decomposition Ro = Ri U Si 
of type (iii) or (iv). Observe that A(S'i) < A(i? ) — ^min by property (3) of 
Lemma [3.61 If e{S\) > eo then procedure for Wo ends; otherwise one iterates by 
applying again EO)l obtaining Si = R2 U S2 with A(S-2) < A(5i) — d m in- It is clear 
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that this procedure cannot continue indefinitely because A(Si) > 0. Therefore, at 
some step e(Si) > eo and the process for Wo ends. 

4. Dependence on the continued fractions 

The computation of ATR points is equivalent to the computation of semi-definite 
integrals of the form 

/TO PC 
/ uj+, ceF. 
J oo 

As we recalled in ij2.21 one can write (|4.ip as a sum of ordinary 4-limit integrals 
by expressing c as a continued fraction with coefficients in Of ■ If the limits of the 
resulting ordinary integrals are too close to the real axis then the number of terms 
to sum for a prescribed accuracy, and therefore the number of Fourier coefficients to 
be computed, may be too large. In this case, the algorithm described in ^Slcan be 
used to express them as sums of integrals whose limits are uniformly bounded away 
from the real axis, reducing the number of terms and Fourier coefficients needed. 

If F is norm-euclidean then the euclidean algorithm gives an effective procedure 
for computing continued fractions. This is the method used in the numerical calcu- 
lations over Q(V29), Q(\/37) and Q(\/41) carried out in }DL03j . But this can only 
be done in a few fields: there are only finitely many norm-euclidean real quadratic 
fields, being Q(\/73) the one having largest discriminant. 

An algorithm for computing continued fractions in 2-stage euclidean real qua- 
dratic fields was given in |GM12j . A field F is said to be 2-stage euclidean if for 
every a, b G Op, b ^ 0, there exist either: 

(i) q,r G Of with a = qb + r and Nm(r) < Nm(6), or 

(ii) Qi, 92, ri, r 2 G Of with a = q x b + n, b = q 2 n + r 2 and Nm(r 2 ) < Nm(6). 

All real quadratic fields of class number 1 are conjectured to be 2-stage euclidean 
(see |Coo76j ). Actually, the algorithm of [GM12j can also be used to verify that 
a given F is 2-stage euclidean, and this was used to prove that all real quadratic 
fields of class number 1 and discriminant up to 8, 000 arc indeed 2-stagc euclidean 
|GM12| Theorem 4.1]. 

Unlike the situation encountered in norm-euclidean fields, 2-stagc division chains 
as in condition ((Ii]) are not unique. As a consequence, elements of F admit in general 
many different continued fraction expansions. This leads to different expressions of 
(|4.ip as a sum of ordinary integrals, whose limits may have very different imaginary 
parts. As it will be illustrated in Sj5]with some explicit examples, numerical exper- 
iments suggest that it is useful to exploit non-uniqueness of continued fractions, by 
searching for continued fractions leading to integrals whose limits have large imag- 
inary parts; or, to be more precise, whose limits give large values of the quantity e 
defined in (|3.2I) . The procedure for computing an expression such as (|4.1[) is then: 

(1) Compute all continued fraction expansions of c given by the algorithm of 
|GM12j . which have length up to a certain fixed bound. 

(2) For each continued fraction, compute the corresponding expression of (|4.1I) as 
a sum of ordinary integrals and compute e m i n : the minimum of the quantities 
e corresponding to the limits. 

(3) Choose the continued fraction giving the highest e m in- 
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(4) For each of the ordinary integrals appearing in the expression given by the 
continued fraction found in the previous step, compute the quantity e and apply 
the algorithm of Theorem 13. II if e < ep, with a suitable choice of e < eo < ep. 

We end the section with some remarks about the algorithm above. 

1. One can exploit non-uniqueness of 2-stage division chains even if F is euclidean. 
As the next section illustrates, this may be beneficial since it usually gives rise 
to integrals with larger values of e, providing an improvement even on the curves 
already considered in |DL03j . In some cases, it may even happen that the e m ; n 
obtained in this way is higher than ep , in which case it is not necessary to apply 
the algorithm provided by Thcorcm l3.ll However, the lack of an a priori estimate 
of the value e m ; n obtained by the non-uniqueness of division chains trick explains 
the key importance of Theorem 13. ll in treating the cases where e m i n < ep. 

2. In Step §3§ we choose the continued fraction giving the highest e m i n because 
experimentally this seems to produce the fewer resulting integrals in step (4). 
We have no rigorous explanation for this fact, although it seems reasonable that 
better initial conditions give better results. 

3. There is a trade-off between small and large values of eo in Step (j4|) above: 
smaller values yield less integrals after the breaking process, but each of these 
integrals requires more Fourier coefficients at the time of integration; on the 
other hand, higher values lead to integrals requiring less Fourier coefficients, but 
the number of resulting integrals tends to be higher. Experimentally, wc found 
that the running time of the algorithm is more sensible to the number of needed 
Fourier coefficients, so a value of eo close to ep seems to be a good choice. For 
instance, in the implementation of the algorithm used to compute the numerical 
examples of the next section, we used eo = 0.81eF, which corresponds to a value 
of ei = 0.9ep. 



5. Numerical verification of Darmon's conjecture 

In this section we illustrate the algorithm described above by calculating approx- 
imations to ATR points which add numerical evidence on top of the one presented 
in [DL03] . Before detailing the computation of an ATR point on -E509 , we comment 
on some calculations of ATR points on three Q-curves that we denote E29, -E37, and 
Eiog . The curves -E29 and -E37 were also considered in |DL03| , and we have included 
them here in order to compare the computational requirements of the algorithm 
used in [DL03] with the one presented in this note. The curve -E109 is an example of 
a curve of conductor 1 defined over a real quadratic field of class number 1 which is 
not norm-euclidean. Therefore, it was not numerically accessible before, although 
it is a Q-curve and algebraic points can be more efficiently computed by using the 
Heegner point method of [DRZ12] . 

The computations for -E29, -E37, and -E109 were performed on a laptop with Intel 
Core™ i5-2540M CPU running at 2.60 GHz, and 8 GB of memory. For the curve 
-E509 w e used a machine equipped with eight Quad-Core AMD Opteron™ Processor 
8384 f° r a total of 32 cores running each at 800 MHz, and equipped with 320 GB 
of memory. 
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The curve E29. Consider the curve defined over F = Q(v29) and given by the 
equation 

1 + v / 29 

E 29 : y 2 + xy + {5uj + ll)y = x 3 , w = . 

For this field the estimated Cf is Cf — 5 (see Remark |3.3[) , which by (|3.6[) yields 
e F ~ 0.0736. 

We consider the ATR field K = F(0) with (3 = y/9uj + 3, for which E 29 {K) 
has a non-torsion point with x-coordinate equal to —1/3. With the algorithm used 
in [DL03] . one obtains integrals with e m i n ~ 0.00145. In order to get 12 decimal 
digits, which is the minimum precision in which the calculations in [DL03 were 
performed, one would have needed to find the Fourier coefficients of all ideals up 
to norm N ~ 6.7 • 10 7 . 

Using the non-uniqueness of continued fraction expansions as explained in Sec- 
tion 31 considering expansions of length up to 5, we obtained 5 integrals with 
imaginary part e m j n ~ 0.0072, which is almost 5 times better than before. In order 
to obtain the same precision one would have to find the Fourier coefficients of ideals 
up to norm TV ~ 2.7 • 10 6 , which is almost 25 times less. 

Since e min < ep, we broke further the integrals with the algorithm of Thcorcm l3.ll 
to move the imaginary parts of the limits close to this theoretical optimal, with a 
choice of eo = 0.81ep in Step (4) of the algorithm outlined in Sectional This yielded 
539 integrals with an imaginary part of e m i n ~ 0.0596, and allowed us to obtain 
the same precision of 12 digits by only considering ideals up to norm TV ~ 40, 000: 
an improvement by a factor of 1,675. By taking ideals of norm up to 40,000 we 
obtained that 

J T = 13.2923360157968468468 ... - 10.78402031269077180934 ...i, 

and —3 • J T coincides with P up to the prescribed accuracy. The calculation took 
less than two minutes. 

The curve E 37 . Let E 37 be the curve defined over F — Q(\/37) appearing in |DL03j 
and having equation 

1 ~\~ \/37 

E 37 : y 2 + y = x 3 + 2x 2 - (19 + 8lo)x + 28 + llw, u = ^ — . 

For this field the constant Cf is approximately equal to 6. By p.6[) we see that 
e F ^ 0.044. 

We consider the field K = F(/3), with f3 = \JAlq + 10, and one of the points 
computed by [DL03] . namely 

P = {-j3 2 /8 - 3/4, -/3 3 /8 - 1/2). 

Using the algorithm of |DL03j one obtains a minimal imaginary part of e m in — 
0.0012, which means that one has to integrate using the Fourier coefficients up to 
norm N ~ 1.12 • 10 s for obtaining 12 digits of precision. 

In order to illustrate the algorithm of Theorem 13.11 we rewrote the integrals 
provided by the method of |DL03] as a sum involving 328 integrals (with a choice 
of e = 0.81cf). The minimal imaginary part then improved to e m ; n ~ e ~ 0.0359, 
which means that to get 12 digits of precision it is enough to use ideals of norm up 
to TV ~ 138, 000. This is an improvement by a factor of 815. In this case, it took 
about 7 minutes to find that 

J T = - 1.3589031642485772101... + 8. 36575277665729384437... i, 
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which satisfies the equality 

5J T = -8P 

up to the prescribed accuracy. 

The curve £109 • In the two remaining subsections we present larger examples 
that were not available to [DL03] . First, consider the curve £109 defined over the 
field F = Q(V109), and given by the equation 



£109 : y 2 + uxy = x 3 — (1 + uj)x 2 - (58cj + 245)x - 630cj - 2944, 

Although £109 is a Q-curve, the field £ is not norm-euclidean and therefore this 
example was not available before. For this field we have that Cp — 10.4, giving that 
ep ~ 0.006. By using the algorithm of Theorem 13. II with, say, eo = 0.81eF, one can 
express any integral U)~f as a sum of integrals with e ~ 0.81 • ep ~ 0.0048. 

In order to compute any such integral with e = 0.0048 to a precision of 12 digits, 
one needs to sum the Fourier coefficients of norm up to roughly 2 • 10 7 . 

Let us consider the point P = (3w + 11, i/3 — 7w — 81/2) defined over the field 
K = £(/?), with f3 = v / 26&ir+T265. Exploiting the non- uniqueness of continued 
fractions, and considering expansions of length up to 5, we obtained 8 integrals 
with e m j n ~ 0.035, which is roughly 6 times higher than ep. It is not necessary then 
in this case to further break the integrals. We computed an approximation to the 
ATR point by considering ideals of norm up to TV = 430, 000, obtaining that 

J r = -3.24024368505944150 . . . • 10" 12 - 42.392087963225793791 ...i, 

which satisfies 

J T = -IP 

up to the prescribed precision of 12 digits. This computation took less than 3 
minutes. 

The curve £509- We consider here the curve already mentioned in Section 1 and 
in |DL03j . defined over £ = Q(v / 509) and given by the equation 

y 2 - xy - ujy = x 3 + (2 + 2w)x 2 + (162 + 3w)x + (71 + 34w), oj = 1 + v ^ . 

We have that Cf ~ 22.5, and therefore tp ~ 0.0015. Theorem 13.11 allows us to 
express any integral Jj° JJ° ojJ as a sum of integrals having e < e < tp . For 
instance, for a choice of eo = 0.81ep we obtain that each of those integrals could 
be computed to 12 digits of precision by summing over the Fourier coefficients of 
norm up to roughly 1.6 • 10 9 ; a bound which, although being large, is within reach 
of the current technology. 

We consider the ATR field K = F{J3) where (3 = ^9144^ + 98577, and the point 
P € E(K) of infinite order having coordinates 

P = {uj + 17, P/2 + V509/2 + 9). 

The field extension A'/£ has relative discriminant of norm 55, which is relatively 
small. Write 

O k =Op + aO F , a 2 + a = 127v / 509 + 2865. 



>K = -t- aUF, 

The ATR points are conjectured to be defined over the Hilbert class field of K. 
Since K has class number 2, we will need to compute the points corresponding to 
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the two non-equivalent optimal embeddings if we want to obtain a point defined 
over K. The first of these embeddings maps 

'0 254a; + 211 
■1 



, , fO 254a; + 2738\ 



whereas the second maps 



a i ^ 952(a) = 



f 27a; + f 369 



The fixed points for the induced action of K > 
vq : K C are, respectively: 



1 



on 'H given by the embedding 



J2) 



0.5 + 0.024492046328012136937583 . 



1 



Exploiting the non-uniqueness of quadratic continued fractions we obtain 4 integrals 
for the first of the points, and 8 for the second. The then minimal imaginary parts 
are e^ n = 0.01917 and e^ n = 0.002926. Since > e F and ej^ n > e F we see that 
in this case it is not necessary to break further the integrals using Theorem 13. II 

In order to obtain about 12 decimal digits of accuracy we precomputed the 
Fourier coefficients of all ideals up to norm 4 • 10 s . The total computation time 
was under two days on the 32-processor machine specified at the beginning of this 
section. We should note that in this computation we heavily exploited parallelism, 
both when computing the Fourier coefficients as well as during the integration step. 

are 



The period lattices for Ek attached to the Neron differential u>e k 
A = (-5.38425378853615683456 . .., -7.44383552310672504690 . 
Ai = (2.47855898378449003059..., 1. 14589256545011559322... i) 
A preimage of P on C/A under the Weierstrass map is 



dx 



2y—x—tu 
A ) = ( A J' A o' 

= (a+ Ar). 



z = -2.6921268942680784172834- • 
The computed values are 
4 1 ' = 22.63291528772669504213102498 . 
A 2) = 106.761524788388057098773188 . . 



5.1426086531573572370822 ...i. 



.1, 



Setting J T 



Jt 2 ^ we find that: 



J - A 



10A^ 



73.6179507973347981534240251 



5.126285-10" 



which suggests that 



J T = Az. 
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